Eigenfrequencies of two mutually interacting gas bubbles in an acoustic field 



O . 

o ■ 

(N ; 

5-H ! 
Oh- 

<\ 

in : 

V • 
i : 

^: 

C/3 ■ 

O ■ 

• i-H . 

C/3 . 
>v 



(N 
> 

O 

oo 
o 



C/3 

o : 
43 : 

Oh. 



X 



Masato Ida 

Satellite Venture Business Laboratory, Gunraa University, 1-5-1 Tenjin-cho, Kiryu-shi, Gunma 376-8515, Japan 

E-mail : ida@vbl.gunma-u.ac.jp 
(February 2, 2008) 

Eigenfrequencies of two mutually interacting gas bubbles in an acoustic field are discussed theo- 
retically and numerically. It is shown by a linear theory that a bubble interacting with a neighboring 
bubble has three eigenfrequencies that change with the distance between two bubbles, and the sign 
and magnitude of the primary Bjerknes force acting on bubbles also change according to the change 
in the eigenfrequencies. 
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INTRODUCTION 



It is known that many bubbles levitated in a liquid 
form a stable structure named "bubble grapes" when a 
weak standing sound wave is applied jl] [3J . In this stable 
structure, bubbles or clusters of them maintain a distance 
each other and tend not to collide, contrary to predictions 
given by the classical theory by Bjerknes which allows ei- 
ther only attraction or repulsion. The distance between 
bubbles or clusters in the structure is comparable to their 
sizes. A similar phenomenon can be observed in a strong 
acoustic field || in which bubbles form a complex fila- 
mentary structure varying slowly, whose scale is much 
smaller than the wavelength of the acoustic field. 

In order to understand these phenomena, one should 
know thoroughly how the mutual interaction of the bub- 
bles influences their acoustic properties. In Ref. |Q, 
Zabolotskaya showed theoretically that the effective reso- 
nance frequencies of two interacting bubbles rise or fall as 
bubbles approach each other and that the variation in the 
resonance frequencies may sometimes causes the change 
of the sign of the secondary Bjerknes force, which is an 
interaction force acting between pulsating bubbles, at a 
certain distance. In Ref. 0, Doinikov and Zavtrak intro- 
duced a similar result by taking the multiple scattering of 
sound between two bubbles into account. They explain 
the origin of the change in the sign of the force as fol- 
lows: When two bubbles oscillate in phase, their surfaces 
move against each other resulting in a stiffening of the 
bubbles. This may increase the effective resonance fre- 
quencies (or eigenfrequencies) of both bubbles and may 
sometimes cause the change in the sign of the force from 
attractive to repulsive. 

Influences of mutual interaction have also been studied 
in regard to sonoluminescence ||,fj]]. In Ref. ||, Mettin 
et al. examined numerically the sign and magnitude of 
the secondary Bjerknes force acting between two small 
air bubbles in a strong acoustic field. They determined 
that even when the distance between two bubbles is much 
larger than their equilibrium sizes the magnitude of the 
force becomes much larger than that predicted by a linear 
theory, and the sign of the force changes near the dynam- 



ical Brake threshold or around the nonlinear resonance 
frequency. In Ref. ]7j], Doinikov showed that a neigh- 
boring bubble drastically varies the magnitude and sign 
of the primary Bjerknes force, which an incident sound 
wave exerts on a bubble, in a strong acoustic field. 

In this paper, following Zabolotskaya's approach, the 
influences of the mutual interaction on the effective res- 
onance frequencies of bubbles are theoretically examined 
in detail. As was pointed out and was mentioned pre- 
viously, the variations in eigenfrequencies strongly affect 
the primary and secondary Bjerknes forces. This means 
that full understanding of the eigenfrequencies surely 
must be of help in knowing how the forces are changed by 
mutual interaction. In Sec. II A, we show that the varia- 



tion in the eigenfrequencies with respect to the distance 
between two bubbles is more complicated than those ex- 
pected by Za bolo tskaya Q and Doinikov and Zavtrak 
P,p|. In Sec. II B, additionally, we discuss theoretically 



the influence of the mutual interaction on the amplitude 
of the bubble pulsation. In Sec. Ill, we show some nu- 
merical results of the primary Bjerknes force, acting on 
two interacting bubbles, by using a nonlinear model in 
order to verify the linear theory for the eigenfrequencies. 
When the frequency of an external sound is equivalent to 
an eigenfrequency of a bubble, the force vanishes because 
of the phase difference of ir/2 between the sound and the 
bubble pulsation. Namely, by observing the points where 
the force vanishes, one can know how the eigenfrequen- 
cies changed by the interaction. 

In the present study, an external sound field of weak 
or moderate amplitude with a frequency comparable to 
the partial resonance frequencies of bubbles is assumed. 



I. MATHEMATICAL FORMULATION 

A gas bubble levitated in a liquid pulsates when a 
sound wave is applied. The sound pressure at the bub- 
ble position drives the pulsation. When other bubbles 
(named "bubble 2" ~ "bubble N" , where N is the num- 
ber of existing bubbles) exist near the bubble (bubble 
1), the pulsation of bubble I is also driven by the sound 



1 



wave scattered by the other bubbles. Namely, the driving 
pressure acting on bubble 1, p$ i , is expressed as 



N 



Pdi = pcx + y ,Psij 



J=2 



(1) 



where p G x and p s \j are the sound pressures of the external 
sound field and the scattered wave emitted by bubble j, 
respectively, at the position of bubble 1. The pressure of 
the external sound can be considered uniform when the 
wavelength of the standing wave is much larger than the 
radii of the bubbles and the distance between bubbles. 
By assuming that those bubbles keep being spherically 
symmetric and the surrounding liquid is incompressible, 
the scattered pressure can be estimated with |6| 



Pslj 



P 



(2) 



rij dt 

where p is the density of the liquid, rij is the distance 
between the centers of bubble 1 and bubble j, Rj is the 
radius of bubble j, and the dot denotes the time deriva- 
tive. 

When the compressibility of the material surround- 
ing the bubbles is completely negligible, one can use the 
RPNNP equation [§: 
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with 
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for modeling the radial pulsation of bubbles, where Pq is 
the static pressure, a is the surface tension at the bubble 
surface, p is the viscosity in the liquid, k is the poly- 
tropic exponent of the gas inside the bubbles, and R\q 
is the equilibrium radius of bubble 1. Substitution of 
Eqs. (0) and (§) into Eq. (|) yields 
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When the number of bubbles is two, this equation is re- 
duced to 
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(4) 



where D = r\i (= r%i). Exchanging 1 and 2 in the sub- 
scripts yields the model equation for bubble 2: 
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with 
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In the following sections, we present theoretical and 
numerical discussions regarding the above systems of 
equations ^ and (||). 



II. LINEAR ANALYSIS 



In this section, some basic properties of Eqs. (Q) and 
(||) are discussed by means of linear theory. By as- 
suming that the bubble pulsation can be represented 
as i?i = R w + e_u R 2 = R20 + e-2, and \ei\ <C Rw, 
I e 2 1 i?20, Eqs. (Q) and ((H) are reduced to the following 
linear formulae: 



ei + w 10 ei + <5iei = 



&2 + w 20 e 2 + <^2e2 
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where 



win = , / — ^j-[3kP + (3k - 1) ^-] 



pRio 



Rw 



and 



^20 



1[3 K P + (3/^-1)!^] 
pn 2Q K20 



are the partial resonance (angular) frequencies of bubble 
1 and bubble 2, respectively, 
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In the case where a ~ 0, the system of equations (p) and 
(0) corresponds to that derived by Zabolotskaya on 
the basis of Lagrangian formalism. 



A. Eigenfrequencies of the linear system of equations 



Let us analyze in detail the eigenfrequencies of the lin- 
ear coupled system, although a brief discussion is shown 
111 Ref. f| . When only one bubble exists in an acoustic 
field and the frequency of the external sound is equiva- 
lent to the resonance frequency of the bubble, the phase 
difference between the bubble pulsation and the exter- 
nal sound becomes 7r/2, as is well known. According 
to this analogy, we call the frequency of an external 
sound an eigenfrequency of the linear coupled system 
when the phase difference between p ex and e± (or p ex and 
e^) becomes ir/2. We assume that the external sound 
pressure at the bubble position is written in a form of 
p cx = —PaSinut. A harmonic steady-state solution of 
the linear coupled system is 
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ei = Ki sm(wi + 4>x), 



(8) 



ffiF = 0, 



where 



or 
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Here the solution for only bubble 1 is shown. The phase 
difference </>i becomes it/2 when 



A x = 0. 



(12) 



It should be noted that a case in which both A\ and B\ 
become zero does not exist. From Eqs. ( |l0| ) and (11), 
one obtains 



H\ 



A " + B " ~ F 2 + G 2 



-RlO/O 



The numerator of this equation always has a nonzero 
value since M 2 > 0; this result denies the existence of 
the case where both Ax — and Bx = are true. It 
should be noted also that F 2 + G 2 appearing in the dom- 
inator of Eq. ( |To| ) always has a nonzero value. When 
G = 0. F is reduced to 
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This has a nonzero, negative value because of 
M 2 L 2 /Mx > 0, RxqR^/D 2 > and MiM 2 > 0. This 
result means that no case exists where both F = and 
G = are true. In a consequence, Eq. (|lj) is reduced to 



HxF + M 2 G = 0. 



(13) 



In the following, we analyze this equation. 

When the viscous terms in Eq. (|1J) are negligible (but 
exist), one can easily obtain the solutions of this equa- 
tion. By assuming that Mi as and M 2 as 0, one obtains 



and 



Hx — L 2 



R 



20 2 



0. 



(15) 



Equation (Q) corresponds to that given in Ref. ||, and 
is symmetric; namely, exchanging 1 and 2 (or 10 and 20) 
on the subscripts of variables in this equation yields the 
same equation. This means that two bubbles have the 
same eigenfrequencies. 

The solutions of Eq. (Ill) are 
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Those solutions can be rewritten approximately as fol- 
lows in the case where the distance D is large enough 
and uixo > ^20'- 
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(17) 
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For D ~ 00, the former and the latter one converge 
to uj 2 q and lo 2 q, respectively. As the bubbles approach 
each other, the former increases and the latter decreases; 
namely, the higher eigenfrequency becomes higher and 
the lower one becomes lower. 

In the case of win < uj 2 o and large D, the solutions 
(HQ) are reduced to 



and 



J 20 
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(20) 



Same as the result for u>xo > ^20? the higher eigenfre- 
quency becomes higher and the lower one becomes lower 
as D decreases. 



The solution of Eq. (15) given for the first time is 

, )2 

,2 _ W 20 



1-R 20 /D' 



(21) 
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This converges to w|o f° r an infinite D and increases with 
decreasing D. In contrast with Eq. (|l4|), Eq. (|l5| ) is asym- 
metric; namely this serves to break the symmetry of the 
Zabolotskaya's result mentioned above. 

The present results show that, when the viscous effect 
of the surrounding liquid is negligible and the radii of the 
bubbles are not equivalent, the bubbles have three asym- 
metric eigenfrequencies; one is the "fundamental eigen- 
frequency" that converges to the partial resonance fre- 
quency of a corresponding bubble for an infinite D, and 
the remaining two are the "subeigenfrequencies" that 
converge to the resonance frequency of a neighboring 
bubble. One of the subeigenfrequencies always increases 
as bubbles approach each other. The other subeigenfre- 
quency decreases (increases) and the fundamental eigen- 
frequency increases (decreases) when the partial reso- 
nance frequency of the bubble is higher (lower) than that 
of the neighboring bubble. 

Now we discuss briefly the case of identical bubbles. 
When two bubbles have the same radii and pulsate in 
phase, Eq. (||) is reduced to 

1 , ^10 \ •• , 2 , x ■ -Pcx 

1 + — jex+^ex+fxe^- — . 
This equation has only one eigenfrequency j^j of 



this difference becomes stronger as D decreases. How- 
ever, all models predict the same result in a qualitative 
sense, which reveals the downward shift of the eigenfre- 
quency. 




1 + Rw/D* 



(22) 



which converges to uj\q for infinite D and decreases with 
decreasing D. Using this example, we here carry out 
quantitative comparison between the model used here 
and previous ones. In Ref. |Io| ], Kobelev and Ostrovskii 
discussed the eigenfrequency of two interacting bubbles 
by means of a mathematical model in which the sur- 
face tension and viscosity are neglected and the poten- 
tial field forming between the bubbles is accurately de- 
scribed by the method of images. The model predicts 
only two eigenfrequencies when Rio ^ R20, as in Zabolot- 
skaya's model, and gives u>\ — Vm 2 uiio ~ 0.832 uj\q 
when Rio = R20 and both bubbles make contact with 
each other, i.e., D = 2i?io- In Ref. |ll]], Zavtrak also dis- 
cussed the eigenfrequency of the present case by means 
of a mathematical model in which multiple scattering 
between two bubbles is taken into account by the Leg- 
endre polynomial expansion. In the case of D = 2i?io, 
Zavtrak's model gives uii ~ 0.832 u>io, which is in close 
agreement with the result obtained by Kobelev and Os- 
trovskii. However, Eq. ( Bp ) gives cui ~ 0.816 uiiq which 
is lower by about 2 % than these values. This quanti- 
tative discrepancy is caused from the difference of func- 
tions for describing the scattered sound wave by a bub- 
ble. The amplitude of the scattered wave is assumed to 
be inversely proportional to D in the model used here 
(See Eq. (§)), while those in the models in Refs. @jn) 
are approximated by an expansion function. Namely, the 
amplitude of the scattered wave at the position of a neigh- 
boring bubble is estimated differently. The influence of 



FIG. 1. The eigenfrequencies of bubbles 1 (cui) and 2 {u>2) 
for Ro ~ 10 /im and fi ~ normalized by lo\o- The dashed 
line denotes the line of ujj /ljio = 1. 
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FIG. 2. The eigenfrequencies for Ro ~ 10 /jm normalized 
by Uio. The dashed line denotes the line of Wj/iOio = 1- 

Figure [l] shows the eigenfrequencies as a function of 
I = D/(Rio + R20), in the case where the viscous effect 
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is negligible. Here, the parameters are set to p = 1000 
kg/m 3 , P = 1 atm, a = 0.0728 N/m, \i w 0, and k = 1.4. 
Equations (M) an d (21j ) are used in plotting those graphs 
since Eqs. (|17j)— (|20) provide a less-accurate result for 
small D. The radius of bubble 1 is fixed to 10 /im and 
that of bubble 2 varies from 5 /im to 15 /im. The eigen- 
frequencies shown in those graphs are normalized with 
the partial resonance frequency of bubble 1, lj%o- As dis- 
cussed above, three eigenfrequencies changing with the 
distance between bubbles appear in each graph except for 
the case of R w = R20 where only one decreasing eigen- 
frequency appears. It should be noted that the larger 
bubble, whose partial resonance frequency is lower than 
that of the smaller bubble, has the highest eigenfrequency 
among those of a pair of bubbles. This result is inconsis- 
tent with the opinion represented in Refs. |^,|| in which it 
is assumed that the effective resonance frequency of the 
smaller bubble rises faster than that of the larger bubble 
when two bubbles approach each other. 

Now we present numerical solutions of Eq. (|l^) for 
examining the influences of viscosity on eigenfrequency. 
Figures |^— show the results obtained by using /i = 
1.137 x 10 kg/(m s) that corresponds to the viscos- 
ity of water at room temperature. From those figures, 
we can observe that, as the viscous effect grows strong, 
i.e., the mean radii of the two bubbles become smaller, 
the subeigenfrequencies vanish gradually from the large- 
distance region (and sometimes from the small-distance 
region) and only the fundamental eigenfrequency re- 
mains. In the case of Ri = 1 /im, the subeigenfrequencies 
of a larger bubble disappear, and when the sizes of two 
bubbles are close each other the fundamental eigenfre- 
quency and the higher one of subeigenfrequencies of a 
smaller bubble vanish in the small-distance region. 

In the case of R\ = 0.1 /im, it is difficult to distinguish 
the subeigenfrequencies from the fundamental one since 
only a smooth curve, which decreases monotonically with 
decreasing I, appears in the graphs. (The eigenfrequency 
of the smaller bubble remaining in the small-distance re- 
gion may be the lower one of the subeigenfrequencies.) 
This result means that in a high-viscosity case the ef- 
fective resonance frequencies of both bubbles decrease 
monotonically as / decreases. 

These results show that the influence of the mutual in- 
teraction on the eigenfrequencies depends on the viscos- 
ity of the surrounding material and is weakened as the 
viscous effect becomes stronger, i.e., the threshold of the 
distance for the appearance of the subeigenfrequencies is 
shortened. 

From the results shown in Fig. ^, we can expect that 
when i?i = 10 /im, Ri — 8 /im, and u> — 1.4cjio, for 
example, the signs of the primary Bjerknes force acting 
on bubbles 1 and 2 change twice as the two bubbles ap- 
proach each other since the curves shown in the figure for 
this case close with the line ujj/ujw = 1.4 (j = 1, 2) at 
two points. We can expect also from the results shown 
in Fig. |3| that when R\ — \ /im, R2 = 0.8 /im, and 
u) = l.lwio, the sign of the force acting on bubble 1 



changes once during the approach since the curves for 
this case close with the line lui/loiq — 1.1 at one point. 
Such an influence of the mutual interaction on the pri- 
mary Bjerknes force is discussed in Sec. III. 
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FIG. 3. The eigenfrequencies for Ro ~ 1 /im normalized by 
Uiq. The dashed line denotes the line of LJj/uio = 1. 
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FIG. 4. The eigenfrequencies for Ro ~ 0.1 /im normalized 
by W10. The dashed line denotes the line of Wj/wio = 1. 
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B. Amplitude of linear pulsation of two coupled 
bubbles 



Next, we briefly discuss the amplitude of the bubble 
pulsation. By substituting Eqs. and (^fj) into Eq. (|^) 
and rearranging it, one obtains 



Pn 



a = Hi + M| 
1 F 2 + G 2 \R wP 



(23) 



When cj is far apart from the partial resonance frequen- 
cies of both bubbles and the viscous term is negligible, 
Eq. (p3|) is reduced to 



III. NUMERICAL EXPERIMENT 

In this section, we show numerical results obtained by 
using the full system of nonlinear equations [(|J) and (||)] 
and sufficiently small amplitude of sound in order to ver- 
ify the linear theory given in the last section. We in- 
vestigate the variations in eigenfrequencies by observing 
the primary Bjerknes force which is a time-averaged force 
acting on a bubble when the external sound field has a 
spatial gradient. The primary Bjerknes force fill is ex- 
pressed as 



-in 



(Rj v Pox ) , 



(25) 



Kl 



l 2 + r 20 u 2 /d y ( P a 

L t L 2 - R 10 R2 u 4 /Dy \R w p 



Furthermore, neglecting the second order term of Rjq/D 
reduces this to 



L 2 + ggo^Vg 
LiL 2 

R 2 qUJ 2 



Pn 



Rwp 



Pa 

Rwp 



(24) 



As D decreases, the amplitude K\ increases when W20 > 
to, and decreases when u 2 o < lu. Physically, this result 
can be understood as follows: When lu 2 q is sufficiently 
larger than u>, the scattered wave due to bubble 2 os- 
cillates in phase with the external sound. In this case, 
the amplitude of the driving pressure acting on bubble 

1 increases (See Eq. (j|)); thus, the pulsation amplitude 
of bubble 1 is increased. In contrast, when uj 2 q is suffi- 
ciently smaller than u>, the scattered wave due to bubble 

2 oscillates out of phase with the external sound; this 
results in the decrease of the pulsation amplitude. 

Now we consider the case where two bubbles have 
the same radii and the frequency of the external sound 
is equivalent to their partial resonance frequencies, i.e., 
Rio = i?2o and u> = luiq = u> 2 q. This setting reduces 
Eq. (||) to 



1 



{& w <4JD* + 6l)u, 



10 



\RwpJ 



This equation shows that the amplitude of the bubble 
pulsation decreases with decreasing D. This result can be 
explained as follows: The decrease of the eigenfrequency 
represented in Eq. (|2^) causes the difference between the 
sound frequency and the effective resonance frequency; 
resultantly, the amplitude of pulsation decreases. Such a 
decrease in the pulsation amplitude of identical two small 
near-resonance bubbles (Ro = 1.5 /im) in a viscous liquid 
is shown in Ref. Jl^ j by employing a direct-numerical- 
simulation technique. The present theory may explain 
the numerical result. 



where (• • •) denotes the time average over a period of 
the external sound. As was mentioned previously, the 
phase difference between the external sound and the bub- 
ble pulsation becomes ir/2 when the sound frequency is 
equivalent to the resonance frequency of a bubble. In 
this situation, furthermore, the primary Bjerknes force 
vanishes as the result of the inner product by Eq. (|2^). 
This means that we can observe the eigenfrequencies of 
interacting bubbles by examining the conditions under 
which the magnitude of the force becomes zero. 

Here we assume that the external standing wave field 
is in a form of 

Pcx(M) = -P a (z) sin(u>t) 

with P a {z) = P s sin(fcz), where z is the spatial coordinate 
vertical to the pressure nodes, P a (z) is the amplitude of 
the standing wave field at z, and k is the wave number 
of the standing wave field. This assumption yields 



. . cos(A;2;) . . 

\/p cx (Z, t) = k . r ffex(z, t) Z, 
sm b) 



(26) 



where z is the unit vector in the z direction. Substitution 
of Eq. (||) into Eq. (||) results in 



with 



3 sin(fczj) ' 



fj = ( R j Pex(Zj,*)), 



(27) 



where Zj is the position of bubble j in the z coordi- 
nate. The method for computing the force is the same 
as what used in Refs. (6[0]. The full system of nonlin- 
ear equations of (Q) and (||) is solved numerically by 
the 4th-order Runge-Kutta method, and fj is then com- 
puted with Rj given by the numerical computation. The 
time average in Eq. ( [27| ) is done over a sufficiently large 
period after the transition has decayed. Here we as- 
sume that the wavelength of the external sound is much 
larger than the distance between the bubbles, and thus, 

Pex(zi,t) =Pcx(^2,*)- 



G 




FIG. 5. fj ~ I curves for Rio = 10 /im and different values of R 2 o. The sound frequency is set to U) = 0.9ljio (a) and 1.4o»io 
(b), and the sound amplitude is O.OOlPo- The theoretical values of (h, l 2 ), where l\ and l 2 are the normalized distances at 
which fx = and f 2 = 0, respectively, given by Eq. are (a): (1.038, 1.038) for R 20 = 8 ^m, (1.402, 1.402) for R 20 = 9 (im, 
(1.279, None) for _R 20 = 15 and (b): (None, 1.362) for R 20 = 5 /im, (1.703 and 2.330, 1.139 and 1.649) for R 20 = 8 ^m, 
(1.224 and 1.258, 1.089 and 1.166) for R 20 = 9 fj,m. 




FIG. 6. — I curves for Rio = 1 (J,m and different values of R 2 o- The sound frequency is set to uj = 0.95wio (a) and l.lwio 
(b), and the sound amplitude is O.OlPo- The theoretical values of (Zi, l 2 ) given by Eq. (jlj) are (a): (1.574, 1.427) for R 20 = 0.8 
Aim, (2.355, 2.077) for _R 2 o = 0.9 /im, (4.761, None) for R 20 = 1.1 pm, (1.049, None) for R 20 = 1.5 ^m and (b): (None, 2.637) 
for R 20 = 0.5 /im, (None, 2.789) for R 20 = 0.8 /^m, (None, 4.279) for R 20 = 0.9 /im, (1.309, None) for R 20 = 1.5 [im. 
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Figures show fj—l curves for the cases of i?io = 10 
lira., Rio = 1 /im and Rio = 0.1 /zm, respectively. Here 
the system of Eqs. (0) and (0) is used; namely, the com- 
pressibility of the surrounding material is neglected. The 
amplitude of the external sound at the bubble position is 
assumed as P a = 0.001P , 0.01P n and 0.1P for R w = 10 
/im, 1 /xm and 0.1 /im cases, respectively. The frequency 
lj of the external sound is set to 0.9wio and 1.4wio fo r 
Rio = 10 /im, 0.95wio and l.lwio for i?io = 1 /im, and 
0.8wio for Rio = 0.1 /im. In those graphs, we can observe 
some points where the primary Bjerknes force vanishes. 
The normalized distances at the points, given theoret- 
ically by solving Eq. (jl^) as a function of D with the 
given u>, are listed in the captions of those figures. These 
theoretical values are in good agreement with the numeri- 
cal ones represented in those graphs. This result confirms 
our prediction that the sign of the primary Bjerknes force 
changes at distances where one of the eigenfrequencies 
becomes equivalent to the sound frequency. 
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FIG. 7. fj—l curves for Rio = 0.1 fim and different values 
of f?20- The sound frequency is u — 0.8u;io and the sound 
amplitude is O.lPo- The theoretical values of (ii, fa) given by 
Eq. (HH) are (1.279, 1.803) for R 20 = 0.11 /im, (1.474, None) 
for i?20 = 0.15 /im. 



IV. CONCLUSION 

In the present study, the dependency of the eigenfre- 
quencies of two mutually interacting bubbles on the dis- 
tance between the bubbles was theoretically determined 
concretely and in detail. The dependency is more intri- 
cate than those predicted or assumed previously by sev- 



eral authors ErW]. It was shown also that the eigenfre- 
quencies depend also on the viscosity of the bubbles' sur- 
rounding material, unlike the partial natural frequency. 
Moreover, the present theory was verified numerically by 
observing the primary Bjerknes force acting on bubbles. 
These present results may be reasonable for understand- 
ing the reversal of sign of the secondary Bjerknes force 
[§-[|,U and the pattern formations of bubbles or clus- 
ters of the m in an acoustic field jl||§. As was pointed 
out in Sec. II A, Zabolotskaya's theory predicts that both 
bubbles interacting with each other have the same eigen- 
frequencies. This result may give rise to simultaneous 
phase changes of the bubbles, e.g., from in-phase to out- 
of-phase with the external sound, and, as a result, may 
deny the reversal of the sign. The asymmetric eigen- 
frequencies derived in the present paper may be able to 
describe the reversal of the sign more accurately. (This 
subject will be discussed in a future paper.) Further- 
more, these present results will also affect understandings 
on some related subjects such as acoustic localization Q 
and superresonances Jl5[ , which were investigated by em- 
ploying a system containing only identical bubbles. 
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